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Abstract We show that, under certain conditions, the micromaser can act as an ef- 
fective source of highly correlated atoms. It is possible to create an extended robust 
entanglement between two successive, initially unentangled atoms passing through a 
cavity filled with with a nonlinear medium taking into consideration a slight level shift. 
Information is transfered from the cavity to the atoms in order to build up entangle- 
ment. The scheme has an advantage over conventional creation of entanglement if the 
two atoms (qubits) are so far apart that a direct interaction is difficult to achieve. 
The interaction of the atoms with the micromaser occurs under the influence of a two- 
quantum transition process. Interesting phenomena are observed, and an extended 
robust entangled state is obtained for different values of the system parameters. Il- 
lustrative variational calculations are performed to demonstrate the effect within an 
analytically tractable two-qubit model. 
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1 Introduction 

Entanglement usually arises from quantum correlations between separated subsystems 
that cannot be created by local actions on each subsystem. Preparation of quantum 
entanglement between distant parties is an important task required for quantum com- 
munication and quantum information processing [T|. In such processing one usually 
needs to find the entanglement properties and a way to control it, therefore the study 
of the dynamic properties of entanglement is useful for processing quantum informa- 
tion. 

In the last years there has been an intensive research in the field of quantum commu- 
nication that has yielded a variety of methods to distribute bipartite entanglement 
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3,4,5,6,7,8 . Nevertheless, due to the lack of a complete understanding of mixed state 
entanglement and multi-partite entanglement, it is not always clear what is the optimal 
way to distribute entanglement among distant parties. 

A great effort has been devoted to the generation of atomic entanglement and en- 
tanglement between cavity modes through atom-photon interactions 9,10,11,12,13 
14, 15 , 16, 17, 18, 19, 20, 21 , 22 and some notable experimental demonstrations have also 
been performed, for instance Refs. [23l[24ll25l[26ll27ll28ll29ll30| . Of particular interest is 
the generation of entangled states in two-atom systems, since they can represent two 
qubits, the building blocks of the quantum gates that are essential to implement quan- 
tum protocols in quantum information processing. 

A number of studies have shown that entanglement can be created between two objects 
that do not interact directly with each other, but interact with a common field, heat 
bath or thermal cavity field 18,30,31 . The formation of atom-photon entanglement 
and the subsequent generation of correlations between spatially separated atoms have 
been shown using the micromaser 27,29,32,33,34 . 

The micromaser 35,36,37 is appreciated as a practical device for processing informa- 
tion. It stores radiation for times significantly longer than the duration of the interac- 
tion with any single atom [35] . The interaction of an atom with the intracavity field of 
a micromaser will leave the atom-field system in an entangled state. The long cavity 
lifetime implies that the memory of this entanglement can influence the interaction 
with subsequent atoms and nonlocal correlations between these successive atoms can 
be induced leading to a violation of Bells inequality |35| . In other words, the successive 
atoms can interact with the field left by earlier atoms, in this manner gain can be 
seen with the most diffuse of gain media containing, on average, less than one atom at 
any given moment [38; and the dynamics of an atom inside the cavity will modify the 
evolution of the later atoms [34 , 39 . 

For a field that is interacting with more than one atom at a time, the atom-field entan- 
glement was investigated [25] . In this scheme [25] , an atom of fixed position situated 
at a peak of the cavity field becomes entangled with a second atom, at a variable dis- 
tance away from the first atom, via their mutual cavity field with which they interact. 
Furthermore, the entanglement between a pair of atoms pumped at the same time 
through a micromaser has been analyzed in Ref . [27] . In practice it is rather difficult to 
realize such a setup though. The genuine one-atom micromaser, on the other hand, can 
be operated over a reasonably large region of parameter space, and is thus a feasible 
device [34] for generating entanglement between two or more atoms. 
It has become well known that the degree of quantum entanglement depends crucially 
on the physical nature of the interacting objects and the character of their mutual 
coupling. It has been noticed that the above investigations involved mostly the absorp- 
tion or emission of a single photon in an atomic transition. However, involvement of 
more than one photon, in particular, two photons in the transition between two atomic 
levels has been known for a long time [40]. The output radiation from such interactions 
exhibits nonclassical properties such as strong sub-Poissonian photon statistics as well 
as fields with increased photon number [41] , Needless to say, the idea of squeezed light 
has originated from two- photon process [42] , Thus, it would be interesting to study the 
properties of two-atom entanglement in the framework of a two-photon process. More- 
over, the two-photon process introduces a dynamic Stark shift in the atomic transition 
which is related to the magnitude of the electric field of the radiation inside the cavity. 
However, it should be noticed that all these results have been obtained for the case 
where the Kerr medium and the Stark shift are ignored. The Kerr medium 43,44,45 
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can be modelled as an anharmonic oscillator with frequency ui. Physically this model 
may be realized as if the cavity contains two different species of atoms, one of which 
behaves like a two-level atom and the other behaves like an anharmonic oscillator in 
the single-mode field of frequency uj |46| . Such a model is interesting by itself. The 
cavity mode is coupled to the Kerr medium as well as to the two-level atoms. 
A Kerr-like medium can be useful in many respects, such as detection of nonclassical 
states |47] , quantum nondemolition measurement [48] . investigation of quantum fluctu- 
ations [49] . generation of entangled macroscopic quantum states 50151 , and quantum 
information processing [521153] . 

In a previous study [31] it has been shown that the entanglement between two qubits 
that do not interact directly with each other can be created for a very short time after 
the interaction is switched on. 

Our purpose here is to demonstrate analytically the increase of the entanglement time 
created for two atoms (qubits) if the cavity field is filled with a nonlinear medium and 
a slight level shift is taken into consideration. In section II, we introduce our model 
and the obtained wave function that controls the model in a general form at time 
t > 0. A brief discussion of the technique we are going to use to compute the entangle- 
ment (concurrence) of mixed quantum states is introduced in section III. The reduced 
density matrices of some special cases of our final state vector at any time t > de- 
pending on various initial states of the system are computed in section IV, supported 
by discussions of the study. Our conclusions are summarized in section V. 

2 The full system and its solution 

We consider a scheme of the micromaser-type where two two-level atoms traverse a 
high-Q (Q ~ 10 9 ) single-mode cavity one after the other in a manner that their flights 
through the cavity do not overlap 18 ,25 32 ,33 35 ,54 55 ,56 ,57 . There is no direct 
interaction between the two atoms, although secondary correlations develop between 
them. The entanglement of their wave functions with the cavity photons can be used 
to formulate local-realist bounds on the detection probabilities for the two atoms 32 , 
135] . The generation of nonlocal correlations between the two atomic states emerging 
from the cavity can in general be understood using the Horodecki theorem [58]. The 
cavity mode is assumed to be filled with a Kerr-like medium 143,44,45]. Each atom 
has energy levels and \0i) (i= 1,2) such that E\ i — E$ i = Tiujq. We assume that 
the two atoms make individually two-photon transitions of frequency 2ui between the 
nondegenerate states 1 0^ } (the ground state, energy Eq^ i = 1,2) and |1^) (the excited 
state, energy E^). The transitions are mediated by a single intermediate level \k) 
(energy E^, with E\ i > E/. > Eq^; the frequencies for 1 0^ ) — > |fc) and — > \k) are 
lu — A and lu + A, respectively. The frequency ljq includes a spontaneous contribution to 
Stark shift due to a direct dipole transition from the intermediate level \k) to |0i) and 
The coupling constants k\ (for jOi) — > jfc)), ki (for jl^) — > |fe)), and A determine 
the Stark shift parameters /3i and P2 of the two levels and also the coupling k between 
the effective two level atoms, states 1 0^ } and |lj), and the field mode 

Pi = K-lA' 1 , (3 2 = K 2 A~~ , K = K\K<lA~ . (1) 

The atom-field interaction is governed by the Jaynes-Cummings (JC) model via a two- 
quantum process 43I44I59]. It is assumed that the atom-field interaction time is shorter 
than the lifetime of the cavity, so that the cavity relaxation will not be considered. The 
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cavity field is assumed to be filled with a nonlinear medium, namely, Kerr medium 43 
44,45 , while the atoms are assumed to have a shift in their levels due to the interaction 
with the radiation field. Assuming for simplicity the photon mode to be in resonance 
with the atoms, the model Hamiltonian under the rotating wave approximation (RWA) 
reads 

H = uj S 3 + aft A + x A f2 A 2 + - S 3 ) + 8 3 )] 

+«(A t2 S- + A 2 S+); (n = l), (2) 

where ui is the cavity frequency, o^o = 2u is the frequency of two atomic energy level 
difference and k = reiK^Zi -1 is the coupling parameter that connects the field with the 
atomic system. We denote by \ the dispersive part of the third-order susceptibility of 
the Kerr- like medium 43 , 44 , 45 . The operator A' (A) is the field creation (annihilation) 
operator, which satisfies the commutation relation [A, A'] = 1. The operators S+, S- 
and 53 are the usual raising, lowering and inversion operators for the two-level atomic 
system, which satisfy the commutation relations [S3, S±] = ±2S± and [S+, S-] = S3. 
The Hamiltonian given by Eq. (2) can be written in the form 

H= H + Hint 1 (3) 

where Hq represents the unperturbed Hamiltonian that is given by 

Ho = u>0A + S 3 ), (4) 

while the perturbed Hamiltonian is given by: 

//= x ^ t2 4 2 + ii t i[ / 3 1 (l-S 3 )+/3 2 (l + S 3 )] +*(i t2 S-+A 2 S+), (5) 

The state vector \ij)f (t = 0)) of the field is represented by a linear superposition of the 
number state |n), i.e., 

00 

IV>/(* = 0)> = ^F n |n>, (6) 
n=0 

where |n) is an eigenstate of the number operator A^ A = n; A^Ajn) = n\n), and F n 
is, in general, complex and gives the probability of the field to have n photons by the 
relation: 

P(n) = <r#/(t = 0)>(^(t = 0)|n> = l^] 2 . (7) 

As already indicated above, we consider a pair of two-level atoms going through the 
cavity mode one after another. Then the initial state vector of the interacting first 
atom-field system is given by 

00 

\1>a-f(t = 0)) = |Vo(t = 0)} ® = o)> = Fn ^ ® 

where 1 1 x ) represents the state vector of the first atom being in excited state. At any 
instant of time t the joint state vector of the field and the first atom can be obtained 
from the solution of the time-dependent Schrodinger equation 

i^ a -f(t))=H\i>a-f(t)), (9) 
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The time of flight through the cavity t is the same for every atom 18,25,32,33 
35,54,55,56,57 , and the joint state vector of both the two atoms and the field may 
be denoted by \tp a _ a _t(t)), and the corresponding atom-atom-field pure-state density 
operator is: 

p{t) = \rl>(t))W)\; \m) = \^a-a-f(t)). (10) 

In order to quantify the degree of entanglement between the two atoms, the field 
variables must be traced out. One may write the reduced mixed-state density matrix 
of the two atoms after taking the trace over the field variables as: 

Pa-a(t) = Tr field p{t). (11) 

Under the initial condition (6), by solving the Schrodinger equation (9), we obtain 
directly the time-dependent wave function of the atom-field system that evolves ac- 
cording to the form 

oo 

|V»o-/(<)> = F ^ UAn K ll> + Rn+2(t) \n + 2, 0i>], (12) 

n 

with the amplitudes K n (t) and R n (t) given by: 

K n (t)=co S (r n t) + i K [X (2n + l) + n(r \ 1)+2 ^ ]gMM, (13) 



„2 _ 

~ 2r 7', 

and, 



Rn^it) = -iK^/n{n - 1) Sln ^ n - 2t) , (14) 

ln-2 



where T n is given by 



T n = k\I \^(2n + 1) + n(r2 o 1} + ^ 1 2 + (n + l)(n + 2) . (15) 



and A n reads: 



An = K[ X n{n + 1)+ ^ + l)+2r* ^ (16) 
k At 

with r — k\/k2- 

Note that within the delay time between the two atoms the field evolves towards a 
thermal steady state, moreover, repetition of the instant in which the later atoms enter 
the cavity means the same field repeats at this instants precisely when successive atoms 
exit the cavity 30 . 

If an additional atom is prepared in a superposition as 

\i> a {t > 0)) = a|l a > + (a - 1) 2 |0 2 ), (17) 

this atom will interact with the field that has been modified by the passage of the 
first atom. Assuming the flight time t of the two atoms through the cavity to be the 
same, the joint time-evolved wave vector of the tripartite system of the two atoms and 
the cavity system after the second atom leaves the cavity is obtained by solving the 
Schrodinger equation, 

i-|V a -a-/(*)> = H int \i> a _ a _ f (t)), (18) 
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which is expressed as: 
\i>a-a-f(t)) = J2 Fn \ a (e~ 2ltAn [Hn(t)\n,l 1 ,l 2 )+T n+2 (t)\n + 2,1 1 ,0 2 )] 



+ e-^»+ 2 e- <M " [J n + 2 (t)\n + 2, Q u 1 2 ) + K„ +4 (t)|n + 4,0i,0 2 }] 



x (a - l) 2 e -«^ e -^»-= [W„(t)|n, li, 2 > + X n _ 2 (t)\n - 2, l lt 1 2 >] 



+e -2riA„ [y n+2 W | n + 2 , 0l , 2 > + Z„(t)|n, Oi, 1 2 >] J |, (19) 

where the amplitudes ff n (t), 2n+i(t), J n+ i(f), V n+2 (t), W n {t), X n _ 2 (£), F„+ 2 (£), and 
Z n (t) are given by: 

ffn(t) = [^n(t)] 2 , (20) 

r n+2 (i) = ^ n (t)_R n+ 2(t) ; (21) 

J n+2 (t) = K n+2 (t)R n+2 (t), (22) 

V n+i (t) = J R„ +2 Wi? n+4 (t). (23) 

W(n,t) = K n (t)K*_ 2 (t), (24) 

X n _ 2 (t)=tf n (i)_R„(i), (25) 

F n+2 (t) = tf*(*)-R™+2W, (26) 



and 



respectively. 



Z„(t) = [R n+2 (t)] 2 , (27) 



3 Entanglement measure 

For bipartite pure states, the partial (von Neumann) entropy of the reduced density 
matrices can provide a good measure of entanglement. However, for mixed states von 
Neumann entropy fails, because it can not distinguish classical and quantum mechan- 
ical correlations. For mixed states, the entanglement can be measured as the average 
entanglement of its pure-state decompositions Ef(p): 

E f (p) = minY,PiE(i>i), (28) 

i 

with E{ipi) being the entanglement measure for the pure state ipi corresponding to all 
the possible decompositions p — XjjPilV'iKV'il- The existence of an infinite number of 
decompositions makes their minimization over this set difficult. Wooters [16] succeeded 
in deriving an analytical solution to this difficult minimization procedure in terms of 
the eigenvalues of the non-Hermitian operator 



R = pp, 



(29) 
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where the tilde denotes the spin-flip of the quantum state, which is defined as: 

p= (a y ®a y )p*(ay ®a y ), (30) 

where a y is the Pauli matrix, and p* is the complex conjugate of p where both are 
expressed in a fixed basis such as {|e}, \g)}. 

In terms of the eigenvalues of R — pp, Ef(p) (known as the entanglement of formation) 
takes the form 



Ef(p) = H l + l^l-CHp) 
where C(p) is called the concurrence and is defined as: 

C(p) = max ^0, \f\i - \f\^, - \f\~z - \f\^ , 



(31) 



(32) 



with the A's representing the eigenvalues of R — pp in descending order, and, 

H(z) = -slogz- (1- z)log(l - z) (33) 

is the binary entropy. The concurrence is associated with the entanglement of formation 
Ef(p), Eq.(31), but it is by itself a good measure for entanglement. The range of 
concurrence is from to 1. For unentangled atoms C(p) = whereas C(p) = 1 for 
maximally entangled atoms. 

We consider special cases of the initial conditions, namely 

(i) only one atom excited; and (ii) initially both atoms excited. 

We apply two different excitations of the initial field, namely Fock state excitation and 
thermal field excitation. 



4 Special cases 

4.1 case 1. Excition in a Fock state 

If the field is excited in a Fock state, the amplitudes F n in Eq. (6) obey the relation: 

F n = 6 m ,n, (34) 
where m is the photon number of Fock state. 

4-1.1 Only one excited atom 

By setting a — in Eq. (17) and F n = 5 m ,n in Eq. (19), we obtain the wave function 
of the system with field excited in a Fock state and with initially excited atom followed 
by an atom in the ground state. 

Having obtained the wave function of the full system, the corresponding density oper- 
ator of the total system can be easily obtained using Eq.(10). The atom-atom system 
can be described in the basis of product states of the individual atoms 

jli,0 2 ) = |l), (35a) 



|ll,l 2 > = |2), 



(356) 
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|0i,0 2 ) = |3), (35c) 

|0i,l 2 ) = |4). (35d) 

Applying Eq. (11) to obtain the reduced density operator of the two atoms, which can 
be written in this basis as: 

Pa-ait) = PU (t) 1 1) (1 1 + P1A (t) 1 1) <4| + p 22 (t) 1 2) <2| 

+P33(*)|3><3| + P4i(i)|4)(l| + P44(*)|4><4|, (36) 



with 



Pii(t) = \W n (t)\ 2 , (37) 

pu(t) = e —N(2„-i)+4±i] t W n (t)Z*(t) = pUt), (38) 

P22(t) = \X n - 2 (t)\ 2 , (39) 

P33(0 = l^+2W| 2 , (40) 

P4A (t) = l^nWI 2 - (41) 



Expressing the reduced density state (36) in matrix form as: 
Pa-a(t) 



(pit 

P22 
p 33 
\P41 p 44/ 



(42) 



jOa-a(t) 



one may write the spin- flip reduced density state p of p by applying Eq.(30) in the 
form: 

>44 P4i\ 
P33 

P22 k ; 

,P14 pii/ 

From an easy procedure one can obtain the square roots of the eigenvalues of the matrix 
R, given by Eq.(29), which are expressed by the set: 

{V^i} = |\/P22P33, VP22P33, -Re(pi 4 ) + V PllP44 - [^(Pu)] 2 , 

fle(m) - V P11P44 - [Jm(pi 4 )] 2 1 . (44) 

As found above, one may use the largest eigenvalue using Eq.(32) to obtain the con- 
currence C(p) as: 



C(p) = 2 ^y/ P11P44 - [Jm(pi 4 )] 2 - VP22p33^j • 



(45) 



One of the interesting phenomenon described by two-level system is the dynamical 
behavior. In the following, we examine the creation of entanglement in a system that 
consisting of a pair of two-level atoms mediated by quantum field contained in a cavity 
through which the two atoms pass successively. For the case when the initially excited 
atom is followed by the atom in the ground stste, the resulting entanglement, measured 
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Fig. 1 Concurrence C (solid curve) and P22 + P33 (dotted curve) as functions of the scaled 
time Kt. The cavity field start from a Fock state with n = 0.0 where r = 0.0. (a) xl K = 0.0. 
(b) x/« = 0.2. (c) x /k = 1.0 (d)x/« = 2.0 



by the concurrence C, as well as the sum of the populations P22 + P33 ar e depicted in 
figures 1, 2, 3 and 4. We examine the effects of the level shifts as well as of the Kerr- 
like medium on the creation of entanglement between the two atoms mediated by the 
cavity field initially prepared in a Fock state. The case of the effective vacuum is quite 
interesting since in this case C oscillates between zeros and its maximum value (see 
Fig. la). It shows a two-peak periodical behavior with maxima « 0.8 that are reached 
at the maxima of P22 + P33 which also shows periodical behavior with maxima « 0.25. 
In this case , the concurrence C reduces to C = 2^/ pupa where P22 = and P14 is 
always real. In fact, C attains the value zero (corresponding to disentangled atoms) 
when P22 + P33 = (corresponding to atoms in pure state) while strong entanglement 
occurs at P22 + P33 = 0.25 (corresponding to atoms in coherent atomic state). It is 
worth mentioning that a similar behavior was shown in Ref. [54], for one-quantum 
process, but with more oscillations in the same interval of time. As soon as we apply 
the nonlinear medium, an interesting result can be observed. The splitting disentangle- 
ment point of the two peaks begins to disappear gradually by increasing \l K till x/ K 
reaches unity, which shows one peak periodical behavior in approximately same time 
intervals, Figs, lb and lc. Moreover, the same behavior is also observed for P22 + P33- 
In this case, the atomic system exhibits long time intervals of entanglement with strong 
Kerr medium. This is because the fact that the higher values of the Kerr parameter 
allow a complete transmission of the interaction field incident on the atomic system. 
A very strong Kerr medium decreases the entanglement maxima, while the periodical 
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behavior is preserved (see Fig. Id). The above results show that the two atoms exhibit 
long time intervals of entanglement the application of Kerr medium with matching 
Kerr parameter. 

An interesting case is the one when we assume that Ki n\ so that the coupling 
parameter ratio, r = ^ < 10~ 2 , where the effect of one of the coupling parameters is 
very weak (see Fig. 2a). 

We notice that C as well as P22 + P33 evolve identically with fixed-amplitude periodical 
evolution. Moreover, C shows two-peak oscillatory behavior with maxima (~ 0.8) that 
are reached at the maxima of P22 + P33j Fig. 2a. Also, C shows very small rapid oscilla- 
tions around its maxima before it collapses to its minimum (see Fig. 2a). This implies 
a longer time of strong entanglement between the two atoms. When the Stark shift 
parameter r increases, a considerable decrease of the maxima of C is found comparing 
with that of P22 + P33, especially when r = 0.2, similar to the effect of very high Kerr 
parameter, Fig. Id. 

Further interesting results are found when we take into account the effects of both Kerr 
and Stark parameters. One peak periodical behavior with clearly remarkable interval 
( w twice the case when no Kerr and stark) of time with very small rapid oscillations 
around C maxima (Fig. la, 2b). In this case the two atoms are in the entangled state 
for longer periods of time than in the previous cases before collapse to disentangled 
state. When n — 2.0, C falls off rapidly except for some revivals in irregular intervals 
with the minima of C being reached at the maxima of P22 + P33 (Fig. 3a). Moreover, a 
weak stark constant, r = 0.5, reduces the revival intervals with a remarkable reduction 
of the C maxima, Fig. 3b. 

A surprising result is found when the excitation number n = 2.0, while ^ < 10~ 2 . 
The two atoms remain in a pure state forever although the populations sum oscillates 
periodically with time with maxima equal to unity. In this case, due to the strong ten- 
dency of the ground atom to become excited, its interaction with the field dominates, 
while neither interaction between field and the excited atom (since no role of its level 
shifts) nor the two atoms themselves. In this case, the atomic states oscillates between 
li,02} and |li, I2), i- e -> both probabilities pn and P22 contribute to C, while P33, 
P44 and 7m[pi4] are always zero, which implies C = in all times (Fig. 4a). So we 
do not need a full population to get strong entanglement between the two atoms as 
illustrated in figure la which shows the opposite behaviors of P22 + P33 and C, i. e., 
the smaller the the populations sum, the higher the concurrence and stronger entan- 
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scaled time scaled time 

Fig. 3 The same as Fig. 1 but when n = 2.0, (a) \/ n = 0.0. (b) x/« = 0.5. 




Fig. 4 Concurrence C (solid curve) and p22 + P33 (dotted curve) as functions of the scaled 
time Kt. The cavity field start from a Fock state with n = 2.0 where \/ n = 0.0. (a) r = 0.001. 
(b) r = 0.5. (c) X = 2.0, r = 0.1 



glement between the atoms. With increasing r, a kind of entanglement between the 
two atoms is created due to the role of the level shifts of both atoms. It is evident 
that the strongest degree of entanglement occurs when P22 + P33 reaches its minimum. 
Moreover, the amplitudes of both C and P22 + P33 decrease considerably as |r| > 1. A 
case with high Kerr parameter, \ j k — 2.0, with values of r = 0.1 is considered in figure 
4c, which shows the tendency of both atoms to be entangled more quickly with the 



12 



maxima (~ 0.75) being reached when p 22 + P33 = 0.25 and any value of p 22 + P33 l ess 
or greater than this value means low degree of entanglement between the two atoms. 



4-1.2 Two excited atoms 

By setting F n = 5m, n and a — 1, we obtain the wave function of the system with field 
excited initially in a Fock state and with initially the two atoms are excited. 
In this case, the atom-atom system can be described in the basis 



|ll,l2> = |l>, 
|ll,0 2 > = |2>, 
|0i,l 2 ) = |3>, 
|0i,0 2 ) = |4), 



(46a) 

(466) 
(46c) 
(46d) 



where the reduced density operator of the two atoms can be written in this basis as 
Pa-a{t) = Pll (*)|1> <1| +P22(*)|2)(2| + P33(*)|3) <3| 
+P23(*)|2><3| +P32(*)|3) <2| + P 44(*)|4)<4|, 
pu(t) = \H n (t)\ 2 , 
Pzt{t) = \T n+2 (t)\ 2 , 
P23W = e «Pf (2«+3) + ^]t Tn+2 (t)J* +2 (t) = ph(t), 



P33(*) = \Jn+2(t)\ \ 

pa{t) = \V n+i (t)\ 2 . 
One may write p a -a(t) in the form: 



(47) 
(48) 
(49) 

(50) 
(51) 
(52) 



Pa-a(t) = 



/Pll \ 
P22 P23 
P32 P33 
V p 44 / 



(53) 



while the spin- flip reduced density operator p can be obtained by applying Eq.(30) 



Pa-a(t) 



'P44 \ 
P33 P32 
P23 P22 
pii/ 



(54) 



and the square roots of the eigenvalues of the matrix 7?, given by Eq.(29), are the 
following 



{VM = |-Re(p23) + VP22P33 - [^(P23)] 2 ,-Re(p 23 ) - \J p 22 P33 ~ \J' 



m(p 23 )] 2 



VPHP44, VP11P44 



(55) 
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0.2 ; 




Fig. 5 Concurrence C (solid curve) and P22 + P33 (dotted curve) as functions of the scaled 
time Kt. The cavity field start from a Fock state with n = 0.0 where r = 0.0. (a) xl K = 0.0. 
(b) X /k = 1.0. 



By using of Eq.(32) to obtain the largest eigenvalue, the concurrence C(p) takes the 
from 



(56) 



C(p) = max(0, v Ai — y A2 — v A3 — y A4) 



P22P33 - [I™(P23)] 2 - VP11P44 



To obtain a clear understanding of the situation we examine the concurrence and 
population dynamics when both successive atoms enter the cavity initially in an excited 
state. In vacuum, one can clearly notice the strong positive effect of the nonlinear 
medium on the degree of entanglement of the atomic system. When the Kerr parameter 
is absolutely zero, the maximum degree of entanglement (~ 0.75 ) is reached near the 
end of the scaled time (at t ~ 8/n) where the two atoms in some kind of opposite states. 
Moreover, by increasing the value of the Kerr parameter, \/ n = 0.5, the concurrence 
begins with maximum C ~ 0.72 is reached in the near the begin of time scale with 
wider intervals of entanglement of the atomic system (Figs. 5a,b). 
On taking into consideration the Stark shift, the obtained results are illustrated in 
Fig. 6. We can notice clearly the similar behavior as in section (4.1.1) except for 
that C reaches its minimum where P22 + P33 = 1-0 , and when r = 0.2, little small 
shift parameter, a quasi-periodical behavior with wider temporal intervals of atomic 
entanglement are showed due to the Stark shift in opposite to its effect in corresponding 
case. When the cavity is excited with number state n = 2, and \I K — 0-0, xl K = 0-5, 
we notice the strong positive effect of the excitation number on the entanglement of 
the atomic system. We can clearly notice more oscillations of C in the same time 
intervals, see Fig. 7. This implies that the Kerr medium acts as factor of enhancement 
of the entanglement between the two atoms in opposite to the same situation of one 
cavity Fock state |54] , Moreover, the maxima of C depend crucially on the maxima of 
P22 + P33- 

More surprising is the case when ^ < 10~ 2 , in this case, p\\ ~ 1.0, while P22, P33, 
P44 and 7m[p23] are always zero, which implies that the two atoms remain in their 
initial excited states and the cavity field plays no role and entanglement of the two 
atoms is not observed. As the atomic system has a level shift, the atomic system shows 
entanglement whose maxima are reached at the maxima of P22+P33- The entanglement 
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amplitue decreases as the Stark shift parameter increases and as possible as the Kerr 
parameter is still small, see Fig. 8. 

Opppsite to the case of only one excited atom, when n > 0, the number of photons 
in the cavity destroys the entanglement between the two atoms in case of small Stark 
shift. As the level shift between the two atomic levels increases, more intervals of 
entanglement between the atoms are created associated with increasing of the maxima 
of C which occur at the maxima of P22 +P33 till r = 2.0 by which a periodical evolution 
of C appears. 
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Fig. 8 Concurrence C (solid curve) and p22 + P33 (dotted curve) as functions of the scaled 
time Kt. The cavity field start from a Fock state with n = 2.0 where xl K = 0.0. ( a ) r = 0.5. 
(b) r = 2.0 (c) x/k = 0.1, r = 2.0. 



4.2 case 2. Excition in a Thermal state 

The thermal field is the most easily available radiation field. At thermal equilibrium, 
the field has an average photon number given by: 

n={e Hu/kT -I)' 1 , (57) 

with Boltzmann constant k and absolute temperature T. The photon distribution p(ri) 
is given by 

p(n) = — T , (58) 

yy ' (l + n) n + 1 ' v ; 

which has a peak at zero, i.e., n pea fc = 0. 
4-2.1 Only one excited atom 

Setting F n — F n 5 n ,n and a = the wave function that governs the system in a 
thermal state, with initially excited atom followed by the one in the ground stste, can 
be obtained. With the condition that p(n) = \F n \ is the photon distribution function 
of the thermal cavity given by Eq. (58) , the reduced density operator of the atom-field 
system after taking the trace over the field variables has the form of Eq. (36) with the 



16 




scaled time scaled time 



Fig. 9 Concurrence C (solid curve) and P22 + P33 (dotted curve) as functions of the scaled 
time ret. The cavity field start from a thermal state with average photon number n = 0.5 where 
r = 0.0. (a) x/k = 0.0. (b) x/re = 0.5. 



coefficients given by: 

Pll«=^>W|WnM| 2 , < 59 ) 
n 

m(t) = ^pHe-^ 2 ^ 2 "- 1 ^ 1 ^ 1 !* w n (t)z*(t) = pUt), (60) 

n 

P2 2 (t)=$>MI*n-2(i)| 2 , (61) 
n 

P33{t) = ^p{n)\Y n+2 (t)\ 2 , (62) 
n 

P44W = ^p(n)|Z„(t)| 2 . (63) 

n 

With these elements, following the same procedure, one can easily compute the con- 
currence C(p) given by Eq. (32). 

In the following we compare the results obtained when the cavity field is excited in the 
thermal field with various mean photon numbers. The results are depicted in Figs. 9, 
10, 11 and 12. 

A small average photon number, h = 0.5 creates a high degree of entanglement of 
chaotic behavior of the atomic system with many maxima of the highest value (~ 0.88) 
is reached when P22 + P33 = 0.0. Moreover, the atomic system remains entangled for- 
ever, Fig. 9a. A similar effect can be noticed by increasing the average photon number, 
n — 2.0, with higher degree of entanglement C ~ 0.93 accompanied by increase of its 
minima, Fig. 11a. Choosing a Kerr parameter of value x/ K = 0.5, affects the general 
behavior of C negatively, where C goes to zero after one period of P22 + P33, while 
its maxima are remarkably reduced. However, increasing n decreases the maxima of C 
remarkably while the general behavior is preserved, see Fig. lib. 

A Stark shift with parameter r < 10~ 2 creates periodical entanglement with maxima 
(~ 0.82) with period t = 0.7n7r/K, n = 0, 1, 2, ... collapse to minima slower than when 
n = 2.0, see Figs. 10a,12a. 

Increasing the shift parameter, r — 0.1, reduces the maxima of C, while a quasi- 
periodical behavior can be noticed while similarly to the case of n = 2.0, the state of 
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the two atoms is not a pure stste, except for the case when the maxima of C reduce 
remarkably, Figs. 10b, 12b. 

When the effects of both Kerr-like medium and Stark shift are taken into account, the 
only noticeable effect is the reduction of the minima of C while the general behavior 
as the case when no stark is present, Fig. 9c, 12d. 
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Fig. 12 The same as Fig. 10 but for n = 2.0. 



4-2.2 Two excited atoms 

To obtain the wave function of this case, we set F n = F n 8 n ,n and a — 1 in Eq. (19). 
With the same condition p(n) = |-Fn| 2 , the reduced density state of the atom-system af- 
ter taking the trace over the field variables has the form of Eq. (28) with the coefficients 
given by: 

Pll (t) = J2p(n)\H n (t)\ 2 , (64) 

n 

P22(i) = ^p(n)|T„ +2 (t)| 2 , (65) 

n 

P2 3 (t) = ^p(n)e"^( 2 " +3 ' + ^]' T n+2 (t)J* n+2 (t) = ph(t), (66) 

n 

P33(t) = 5>WK+2(i)| 2 , (67) 

n 

P44(«) = ^ P (n) |l/„ + 4 W | 2 . (68) 

n 

With these elements the concurrence C(p) can be easily computed. 
Remarkably interesting results are found when the injected thermal field interacts with 
two atoms passing through it in excited states. The results are shown in figures 13, 
14, 15, and 16. We notice clearly the average photon number reducing the general 
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scaled time scaled time 



Fig. 13 Concurrence C (solid curve) and P22 + P33 (dotted curve) as functions of the scaled 
time Kt. The cavity field start from a thermal state with average photon number n = 0.5 where 
r = 0.0. (a) x/k = 0.0. (b) x/re = 0.5. 




Fig. 14 The same as Fig. 13 but for x/« = 0.0. (a) r = 0.01. (b) r = 0.3 (c) */« = 1.0. 
r = 0.01 (d) x/i* = 0.5, r = 0.3. 



behavior of the concurrence C, while similar behaviors to the corresponding cases of 
Fock state field are noticed. Moreover, the effect of nonlinear medium on increasing the 
maxima of C, and creating a periodical entanglement with small oscillations around 
its maximum with wider periods by taking into account the effect of level shifts are 
preserved. Also, the behavior of C, where reaches its maxima at the maxima of P22+P33 
is also preserved, see Figs. 13-16 and 5-8. 
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scaled time scaled time 



Fig. 15 The same as Fig. 13 but for n = 2.0. 




scaled time scaled time 



Fig. 16 The same as Fig. 14 but for n = 2.0. 



5 Conclusion 

In conclusion, from the results illustrated in the previous sections, we can conclude 
that two atoms (two-qubits) entanglement via two-photon process is more sensitive 
to the initial conditions than one photon process. Long-time intervals of two-qubit 
entanglement can be achieved by filling the vacuum cavity with a Kerr-like medium 
with parameter close to unity and taking into account a slight level shift regardless of 
the initial states of the two atoms. Moreover, when the cavity contains only one atom 
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excited, long-time periods of two-qubits entanglement with no decay to zero can be 
achieved by applying a weak Kerr medium with small level shift when the cavity is 
excited in the thermal state and contains only one excited atom. Furthermore, the two 
atoms become less entangled in excited cavity, while they become stronger entangled 
as well as the effects of both Kerr-like medium and Stark shift taken into consideration. 

Acknowledgements The author would like to thank the referees for their constrictive com- 
ments that helped to improve the text in many ways. 
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